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ABSTRACT 

We present the broad-band 0.6-150keV Suzaku and Swift BAT spectra of the low lu- 
minosity Seyfert galaxy, NGC 7213. The time-averaged continuum emission is well fitted 
by a single power-law of photon index F — 1.75 and from consideration of the Fermi flux 
limit we constrain the high energy cutoff to be 350 keV <£'cut <25 MeV. Line emission from 
both near-neutral iron Ka at 6.39 keV and highly ionised iron, from FeXXV and FeXXVl, 
is strongly detected in the Suzaku spectrum, further confirming the results of previous ob- 
servations with Chandra and XMM-Newton. We find the centroid energies for the emission 
from Fe XXV and Fe XXVI to be 6.60 keV and 6.95 keV respectively, with the latter appear- 
ing to be resolved in the Suzaku spectrum. From modelling, we show that the Fe XXV and 
Fe XXVI emission can result from a highly photo-ionised plasma, with a column density of 
TVh ^ 3 X IQ-^'^ cmT^. A Compton reflection component, e.g., originating from an optically- 
thick accretion disc or a Compton-thick torus, appears either very weak or absent in this AGN, 
subtending < 1 sr to the X-ray source, consistent with previous findings. Indeed the absence 
of Compton reflection from either neutral or ionised material coupled with the lack of any rel- 
ativistic Fe K signatures in the spectrum suggests that an inner, optically-thick accretion disc 
is absent in this source. Instead, the accretion disc could be truncated with the inner regions 
perhaps replaced by a Compton-thin Radiatively Inefficient Accretion Flow (RIAF). Thus, the 
Fe XXV and Fe XXVl emission could both originate in ionised material perhaps at the transi- 
tion region between the hot, inner flow and the cold, truncated accretion disc on the order of 
10"^ — 10"* gravitational radii from the black hole. The origin for the unresolved neutral Fe Ka 
emission is then likely to be further out, perhaps originating in the optical Broad Line Region 
or a Compton-thin pc-scale torus. 
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1 INTRODUCTION 

NGC 7213 is a nearby low-luminosity AGN (z = 0.005839), of- 
ten classified as an intermediate between a Seyfert 1 and a LINER 
(Low-Ionisation Nuclear Emission-line Region) galaxy due to its 
optical spectrum (Filippenko & Halpern 1984). Its X-ray spectral 
properties also appear to lie between those of weak AGN (e.g., 
M81) and 'classical' higher-luminosity broad-line Seyferts. The ul- 
traviolet (UV) flux measured by Wu, Boggess & Gull (1983) was 
higher than would be expected from an extrapolation of the optical 
flux, indicating that NGC 7213 may have a Big Blue Bump (BBB), 
although weak compared to most Seyferts. This object has a high 
black hole mass of about 1O*M0 as estimated from the stellar ve- 
locity dispersion (Nelson & Whittle 1995; Woo & Urry 2002) and a 



low bolometric luminosity (Lboi) of about 9 x lO**^ ergs ^ (Star- 
ling et al. 2005). 

NGC 7213 has a very low accretion rate of only 0.07% I/Edd, 
a value which is intermediate between those usually found in local 
Type 1 Seyfert Galaxies (e.g., Padovani & Rafanelli 1988; Wandel 
1999; Page 2001) and LINERs (e.g.. Ho 1999). Interestingly, this is 
much less than the predicted 2% LEdd "critical" rate whereby the 
high/soft state in black hole X-ray binaries can be observed (Mac- 
carone 2003). Furthermore, this object exhibits another interesting 
characteristic since it is part of a class of Seyfert galaxies which 
have radio properties that are intermediate between those of radio- 
loud and radio-quiet active galaxies (e.g., Blank, Harnett & Jones 
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2005 and reference therein). It is therefore conceivable that NGC 
7213 is an analogue of the Galactic low/hard state sources. 

A simultaneous XMM-Nswton (net EPIC-pn exposure 
~30ks) and BeppoSAX observation in May 2001 revealed fur- 
ther peculiar characteristics of NGC 72 1 3. The low S/N RGS (Re- 
flection Grating Spectrometer) spectrum showed the presence of 
several weak emission features with no absorption lines (Starling 
et al. 2005) contrary to what is usually found in Seyfert Type 1 
Galaxies. Moreover, the emission lines appeared to be the signa- 
ture of a coUisionally-ionised thermal plasma (fceT ^ O.lSkeV), 
while in Seyfert Galaxies, only emission and/or absorption lines 
from a photo-ionised warm absorber / emitter have before been ob- 
served (e.g., NGC 3783, Kaspi et al. 2001; NGC 4151, Schurch & 
Warwick 2002; NGC 1068, Kinkhabwala et al. 2002; Brinkman et 
al. 2002; Mrk 3, Pounds & Page 2005). Such emission lines from 
collisionally-ionised thermal plasma have, however, been observed 
in the soft X-ray spectra of LINERs, such as M81 (Page et al. 2003). 
Interestingly, no significant Compton reflection was observed in 
this simultaneous XMM-Newton and BeppoSAX observation (re- 
flection component: R = Q,/2n < 0.2, where a value R = 1 
corresponds to reflection off material subtending 2-k sr; Bianchi et 
al. 2003), though the presence of a significant Fe K complex could 
be explained by three narrow emission lines: neutral iron at around 
6.40 keV with an EW of ~80eV, Fexxv at around 6.66 keV and 
Fexxvi at about 6.94 keV (see also Starling et al. 2005). 

Bianchi et al. (2003) deduced from the absence of the reflec- 
tion component that the neutral iron Kq emission line is most likely 
produced in a Compton- thin torus or the Broad Line Region (BLR). 
Indeed, according to Matt, Perola & Pirlo (1991) and George & 
Fabian (199 1), a line with an EW of ~80 eV would require a reflec- 
tion component of about R ^ 0.6, a larger value than that found for 
this object {R < 0.2). Furthermore, Bianchi et al. (2008) reported 
the data analysis of a long Chandra/HETG observation of NGC 
7213 finding that the neutral iron Ka line is resolved with a FWHM 
value of 2 400tgoo'' km s" \ fuUy consistent with the Ha hne width 
(2640tgJ°kms"^) measured with the ESO/NTT telescope. They 
therefore inferred that the neutral Fc K lino seen in this object orig- 
inates in the Compton-thin BLR explaining the lack of evidence 
for Compton reflection. They also confirmed the presence of two 
ionised iron lines at ~6.72keV and ~6.99keV which they asso- 
ciate most probably with the resonance transition of the Fexxv 
triplet and the Lya transition of Fe XXVI, respectively. Comparing 
the line energies found with their rest-frame values, a blue-shift of 
about 900kms~^ is inferred. If the dominant line at ~6.72keV is 
indeed the resonance line of the FeXXV triplet, then this means 
that the line may be associated with a collisionally-ionised thermal 
plasma (Porquet & Dubau 2000; Bautista & Kallman 2000). 

Here we report on a 90 ks Suzaku (Mitsuda et al. 2007) obser- 
vation of NGC 7213. The overall goal is to understand accretion in 
an AGN at low rates compared to Eddington, through a high signal- 
to-noise, broad-band observation of this source. Specifically, the 
objectives are to parameterise the iron line complex with an anal- 
ysis of the XIS (X-ray Imaging Spectrometer) spectra in order to 
constrain the individual properties of the lines before considering 
the complete broad-band spectra from 0.6-150 keV with a com- 
bined analysis of the Suzaku XIS and HXD (Hard X-ray Detector) 
data with that from the time-averaged Swift BAT 22-month sur- 
vey (see Section 3). The aims are to obtain better constraints on 
the origin of both the iron line complex and any observed soft ex- 



cess whilst simultaneously testing for the presence (or absence) of 
a Compton reflection component (see Section 4). 



2 SUZAKU ANALYSIS AND DATA REDUCTION 

2.1 Suzaku Analysis 

NGC 7213 was observed by Suzaku on October 22 2006 with a to- 
tal net exposure of 90 ks. In this paper we discuss data taken from 
the 4 XIS (Koyama et al. 2007) CCDs and the PIN diodes of the 
HXD (Takahashi et al. 2007). 

Events files from version 2.0.6.13 of the Suzaku pipeline pro- 
cessing were used. All events files were screened within XSELECT 
to exclude data taken within the SAA (South Atlantic Anomaly) as 
well as excluding data with an Earth elevation angle (ELV) < 5 de- 
grees. Data taken with Earth day-time elevation angles (DYE_ELV) 
less than 20 degrees were also excluded. A cut-off rigidity (COR) 
criteria of > 6 GeV/c for the XIS was applied. Only good events 
with grades 0,2,3,4 and 6 were used, while hot and flickering pix- 
els were removed from the XIS images using the CLEANSIS script. 
Time intervals affected by telemetry saturation were also removed. 

Subsequently, source spectra from the XIS CCDs were ex- 
tracted from circular regions of 2.3' radius centered on the source, 
in the on-axis XIS nominal pointing position. Background spec- 
tra were extracted from 2.3' circles offset from the source re- 
gion, avoiding the calibration sources on the comers of the CCD 
chips. XIS response files (RMFs) and ancillary response files 
(ARFs) were generated using the XISRMFGEN and XISSIMARFGEN 
FTOOLS respectively including correction for the hydrocarbon con- 
tamination on the optical blocking filter (Ishisaki et al. 2007). A net 
XIS source exposure of 90.7 ks was obtained for each of the 4 XIS 
chips. The 3 front-illuminated XIS chips (XIS 0,2,3; hereafter XIS- 
Fl) are predominantly used in this paper as they have the greatest 
sensitivity at iron K. These chips were found to produce consistent 
spectra within the statistical errors, so the spectra and responses 
were combined to maximise signal to noise. The net source count 
rate for the 3 XIS combined was 1.585±0.002 counts s"^ per XIS, 
with the background rate only 0.7% of the source rate. This count 
rate corresponds to an observed flux of 2.46 x 10~^^ ergcm"^ s~^ 
and a luminosity of 1.85 x 10''^ erg s~^ over the 2-lOkeV range. 
The XIS source spectrum was binned at the HWHM resolution 
of the detector due to the high photon statistics. This enabled the 
use of minimisation as there were >50 counts per resolution 
bin. Errors are quoted to 90% confidence for 1 parameter (i.e., 
Ax^ = 2.7) unless otherwise stated. 

2.2 HXD Reduction 

As NGC 7213 is below the detection threshold of the HXD/GSO, 
we used data from the HXD/PIN only, where this object is de- 
tected at the 13(7 level relative to the background. The source 
spectrum was extracted from the cleaned HXD/PIN events files 
and processed with the screening criteria described above. The 
HXD/PIN instrumental background spectrum was generated from 
a "tuned" time dependent model provided by the HXD instrument 
team (Fukazawa et al. 2009). Both the source and background 
spectra were made with identical GTIs (Good Time Intervals) and 
the source exposure was corrected for detector deadtime (which is 
« 6.7%). A detailed description of the PIN detector deadtime is 



© 2010 RAS, MNRAS 000, 1-?? 



The Suzaku View ofNGC 7213 3 



given in Kokubun et al. (2007). The net exposure time of the PIN 
source spectrum was 84.3 ks after deadtime correction. Note that 
the background spectral model was generated with 10 x the actual 
background count rate in order to minimise the photon noise on the 
background; this has been accounted for by increasing the effec- 
tive exposure time of the background spectra by a factor of x 10. 
The HXD/PIN response file dated 2008/01/29 (epoch 3) for the XIS 
nominal position was used for these spectral fits. 

In addition, a spectrum of the cosmic X-ray background 
(CXB) (Boldt 1987; Gruber et al. 1999) was also simulated 
with the HXD/PIN. The form of the CXB was taken as 9.0 x 
10"^(£'/3keV)"''-^^exp(-£;/40keV) ergcm"^ s"^ sr"^keV"\ 
When normalised to the field of view of the HXD/PIN in- 
strument, the effective flux of the CXB component is 8.49 x 
10~^^ ergcm"^ s~^ in the 15-50 keV band corresponding to 
a count rate of ~0.017 counts s~^. The net flux of NGC 
7213 measured by the HXD over the same band is 3.58 x 
10~^^ ergcm~^ s~^, i.e., the CXB component represents 24% of 
the net source flux measured by the HXD/PIN. Note that there may 
be some uncertainty in the absolute flux level of the CXB com- 
ponent measured between missions; for instance, Churazov et al. 
(2007) find the CXB normalisation from INTEGRAL to be about 
10% higher than measured by Gruber et al. (1999) from the HEAO- 
1 data. However, a factor of ±10% uncertainty in the CXB normal- 
isation would result in a ±2.4% uncertainty in the HXD flux for 
NGC 7213, which is well within the statistical uncertainty of the 
HXD/PIN observations. After background subtraction (including 
both the instrumental and CXB components), the resulting net PIN 
source count rate from 15-50 keV was 0.062 ± 0.002 counts s^^ 
corresponding to a 15-50 keV flux of 3.58 x 10^^^ ergcm"^ s~^. 
Note that the total background count rate was ~0.350 counts s~^ 
(15-50 keV) with a typical la systematic uncertainty of ±1.3%. 
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Figure 1. XIS-FI background-subtracted lightcurve of NGC 7213 from 
0.5-10 keV in 5760 s orbital bins. 
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We used 0.6-10 keV data in both the XIS-FI and XIS-BI 
spectra. We ignored the 1.7-1.9 keV band in the co-added FI spec- 
trum and the BI spectrum due to uncertainties in calibration as- 
sociated with the instrumental SiK edge. In all fits, we included 
a constant multiplicative factor to account for relative instrument 
normalisations. We allowed the relative XIS-BI/XIS~FI normali- 
sation to vary, but best-fit values were always within 1% of each 
other. 



Figure 2. Net background-subtracted HXD/PIN lightcurve of NGC 7213 
from 15-50 keV in 5760 s orbital bins. 



sections for X-ray absorption by the interstellar medium from Mor- 
rison & McCammon (1983). Note that all fit parameters are given 
in the rest frame of the galaxy, assuming a distance of 25 Mpc to 
NGC 7213 (Mould et al. 2000). Abundances are those of Anders & 
Grevesse (1989) unless otherwise stated. 



A visual analysis of the lightcurves was undertaken to deter- 
mine whether any detailed timing analysis was necessary. It can 
be seen from Figure 1 that the amplitude of the XIS-FI lightcurve 
varies only by a factor of ~0. 1 throughout the entire observation 
indicating little intrinsic variability below 10 keV. From Figure 2 it 
can be seen that the HXD/PIN lightcurve, too, shows little evidence 
of any substantial variability in the hard X-ray band. Therefore, due 
to the lack of any strong evidence of short-timescale spectral vari- 
ability, we proceed to consider the time-averaged spectrum (Section 
3). 



3 SPECTRAL ANALYSIS 

The XSPEC Vl 1.3 software package (Amaud 1996) was used for 
spectral analysis of the background-subtracted spectrum. In all fits. 



we included the Galactic column density {N^ 



Gal . 



l.I xlO^"cm" 



obtained from the FTOOL NH using the compilations of Dickey 
& Lockman 1990 and Kalberla et al. 2005) and used the cross- 



The cross-normalisation between the HXD/PIN and XIS de- 
tectors was accounted for by the addition of a fixed constant com- 
ponent at a value of 1.16 for the XIS nominal pointing position, 
a value derived using Suzaku observations of the Crab (Ishida, 
Suzuki & Someya 2007^). 



3.1 The Fe Kline profile 

The X-ray spectrum was initially analysed in the 0.6-50 keV band 
using both the XIS-FI and HXD/PIN data. A power-law with 
Galactic absorption of column density A'^e = 1.1 x 10^" cm"'^ was 
fitted to the data revealing a slight soft excess at energies < 1 keV 
as shown in Figure 3. The XIS 1 data were included in this fit to 
check for consistency. For clarity, the HXD data were binned to 
lOcr per spectral bin relative to the background. The hard X-ray 
data are seen to extrapolate quite well to the XIS data with very 



ftp://legacy.gsfc.nasa.gov/suzaku/doc/xrt/suzakumemo-2007-l l.pdf 
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Model Component Fit Parameter Value x^/dof 



1 . Power-law Continuum" 


r 

Normalisa tion 

.F2.5-IO kcV 


1.75 ±0.02 

6 25+°-''2 
— 0.06 

2.18 X 10-11 






2. Galactic Absorption*" 




1.10 X 10^0 


739/316 




3. FeKa Line'^ 


(T 

EW 
FWHM 
Line Flux 


6.39 ±0.01 
< 4.30 X 10-2 
00 1+11.0 

< 4600 

r, 10 + 0. 28 

^•-^"-0.29 


422/313 


316.8 


4. H-like Line'^ 


-E-line 

cr 
EW 
FWHM 
Line Flux 


6.95 ±0.03 

62.3+i|-0 
lOOOO+'^OOO 

iUUUU_4QQQ 


365/310 


57.1 


5. 6.60keVLine= 


cr 

EW 
FWHM 
L ine Flux 


6.60 ±0.03 

< 0.30 
24.4 ±8.0 
< 32000 
0.67 ±0.22 


332/307 


32.8 


6. 8.00keVLine= 


fj 

EW 
FWHM 
Line Flux 


8.00l0;10 

< 0.28 

45.9+27^ 

< 25000 
0.7810.57 


322/304 


10.0 


7. FeK/3 Line= 


(7 

EW 
FWHM 
Line Flux 


7.06 
< 4.30 X 10-2 

< 29.3 

< 4600 

< 0.67 


322/304 


0.5 


8. Fit Statistics"* 


xVdof 
Null Proba bility 


322/304 
0.36 







Table 1. Spectral Pai'ameters to Suzaku XIS in the Fe K band. " r, photon index; Normalisation in units xlO ^photonscm i; F2. 5—10 keV, absorbed 
continuum flux from 2.5-10 keV in units ergscm-2 s-^. ^ Local Galactic absorption (at 2 = 0), units line energy in units keV; a,\a line 

width in units keV; EW, equivalent width in units eV; FWHM, full width at half maximum in units km s— 1; Line Flux in units X 10— ^ photons cm-2 s— 1. 
Reduced chi-squared (x^/rfo/) and null hypothesis probability for spectral fit. 
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Figure 3. The Suzaku spectra of NGC 7213 (in the observed frame) show- 
ing the XIS-FI (black), XIS-BI (red) and HXD/PIN (green). An absorbed 
(Galactic column density) power-law has been fit to the broad-band spec- 
trum. A significant positive residual is observed in the Fe K complex energy 
range, as well as a weak excess in the < 1 keV energy range. 
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Figure 4. A plot of the ratio of the residuals for XIS-FI of the Fe line 
complex compared to the power-law continuum. The spectra are binned 
relative to the half-width at half-maximum of the detector resolution. The 
vertical dotted lines show the expected line energies of, from left to right, 
FeKa, FeXXV forbidden, FeXXV ls-2p resonance, FeXXVI ls-2pLya 
and Fe K/9 in the obsei-ved frame. 
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few residuals in excess of the power-law continuum indicating lit- 
tle or no reflection component {R ~ 0.2, see Section 3.2). Line 
emission is clearly present with a strong but seemingly relatively 
narrow Fe Ka line at ~6.4keV. As the XIS 1 has a much lower 
S/N ratio at higher energies above 2 keV, these data were initially 
excluded from the Fe K line analysis leaving the XIS-FI to be anal- 
ysed from 2.5-10 keV. The HXD data were also initially excluded 
as the lack of Compton reflection suggested that the Fe K emission 
lines could be modelled independently. The HXD and XIS 1 data 
are re-included in the broad-band fits in Section 3.2. 

A simple power-law model of F = 1.75 ± 0.02 with Galactic 
absorption evidently resulted in a poor fit (x^ /dof — 739/316) 
highlighted by a low null hypothesis probability of 6.67 x 10"'^''. 
A plot of the ratio of the residuals with respect to the power-law 
continuum from 5.0-8.5 keV (Figure 4) clearly shows X-ray line 
emission which requires modelling; the most apparent being the 
FeKa line from near-neutral material at 6.39 ± 0.01 keV with 
an intrinsic width of cr < 4.30 x 10~^ keV, equivalent width of 
EW = 83.lt\o° eV and an observed flux of Fkq = ^.IStoH x 
10"^ photons cm~^ s^^. Adding this line improves the fit signifi- 
cantly with a value of Ax^ = 316.8 for 3 parameters of interest. 
However, even upon modelling the strong Fe Ka line at 6.39 keV, 
the fit remains unacceptable (null probability = 3.70 x 10"'') with 
further residuals still present between 6.5-7.0 keV indicating K- 
shell emission from ionised Fe. 

Further Gaussian lines were added to fit other prominent 
emission lines starting with the ls-2p doublet from hydrogen- 
like iron (FeXXVl) at 6.95 ± 0.03 keV which corresponds to a 
value of Ax'^ ~ 57.1 for an additional 3 parameters of inter- 
est. Unlike the 6.39 keV line, this line appears to be resolved 
compared to the detector resolution with an intrinsic width of 
a = O.lOt^EJ^keV (FWHM ~ lOOOOkms"^), an equivalent 
width of EW = 62.3tlH eV and an observed flux of Func = 
1.401q 32 ^ 10~^ photons cm~^ s~^. A third narrow component 
at a line energy of 6.60 ± 0.03 keV was also modelled, improving 
the fit by a factor of Ax*^ = 32.8 for a further 3 parameters of in- 
terest. A line energy vs. line flux contour plot (Figure 5) shows that 
this line energy can be rejected at the 99.9% confidence level (for 
two interesting parameters) as being associated with the resonance 
line of helium-like iron (FeXXV) at 6.700 keV (Ax^ ~ 14) and 
is just acceptable (rejected at only 90% confidence) as the forbid- 
den line at 6.637 keV (also see Section 3.3.2).^ The fact that this 
indicates that the line at 6.60 keV is more consistent with the for- 
bidden transition of Fe XXV is an interesting discovery since the 
resonance line is expected to dominate over the forbidden line in 
a coUisionally-ionised plasma (Porquet & Dubau 2000; Bautista 
& Kallman 2000). Further discussion regarding the origin of this 
emission can be found in Section 4.2. 

A weak, narrow component appeared to remain in the resid- 
uals at a line energy of 8.001q keV with an intrinsic width of 
a < 0.28 keV. However the detection is more marginal; adding 
this extra line component only improved the fit by Ax^ = 10.0 
for 3 parameters. This putative line could possibly be associated 
with the ls-3p transitions of FeXXV or FeXXVI (corresponding 



^ An analysis of the spectrum produced by the XIS Fe-55 calibration 
source, which produces emission lines from Mn Ka and Mn K/3, shows that 
the absolute XIS energy scale is accurate to within ±10 eV. 
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Figure 5. A two dimensional contour plot between line energy and normal- 
isation of the 6.60 keV line at the 68%, 90%, 99% and 99.9% confidence 
intervals (inner to outer contours respectively) for 2 parameters of interest. 
The vertical dashed lines show the rest energies of forbidden, intercombi- 
nation and resonance transitions of the Fe XXV triplet at 6.637, 6.668 and 
6.700 keV respectively. 

to rest energies of ~7.88 and ~8.25 keV respectively) or alterna- 
tively, it could be due to the ls-2p transition of H-like nickel (cor- 
responding to a rest energy of ~8.10keV). Finally, no neutral K/3 
emission was apparent but was still modelled for consistency at a 
fixed line energy of 7.06 keV, with an intrinsic width, cr, tied to 
that of the corresponding Ka line. The upper limit on the flux of 
FKfj < 0.67 X 10"^ photons cm"^ s"'^ then provided an upper 
limit on the Kfi/Ka flux ratio of 0.35. A value of Ax^ = 0.5 re- 
vealed that the the K/3 line is consistent with the fit to the data but 
in this instance is not required. Hereafter we include the K/3 line 
fixed at 13% of the Ka flux in all subsequent fits to remain consis- 
tent with the theoretical flux ratio for neutral iron (Kaastra & Mewe 
1993). 

We did also attempt to model the neutral Fe Ka emission with 
a DISKLINE component (Fabian et al. 1989) to test for the pres- 
ence of any broad, relativistic emission from the inner regions of 
the accretion disc. We fixed the width of the original Gaussian at 
6.39 keV to be narrow (a — 10 eV) in order to model emission 
from distant material and introduced a DISKLINE component to the 
model with the centroid energy fixed at 6.39 keV and the emissivity 
index fixed at a standard value of g = 3. We also fixed the inner 
and outer radii of the emission at 6 and 400 Rg from the black hole 
respectively (where 6Rg corresponds to the innermost stable or- 
bit for a Schwarzschild black hole) and the inclination angle of the 
source at ^ = 30° . Upon fitting, this returned a value for the flux of 
the line of -Fdiskiinc < 7.64 x 10~^ photons cm~^ s~^ which cor- 
responds to a 90% upper limit on the equivalent width of the line of 
EW <26.8 eV. This tight constraint appears to exclude the presence 
of any Fe K emission from the inner accretion disc. We note that no 
other significant emission or absorption lines are found in the Fe K 
band. The values of all of the final parameters and fit statistics are 
noted in Table 1. 

3.2 The Broad-Band Spectrum 

The first stage of the broad-band spectral analysis was to model 
the spectra above lOkeV by including the Suzaku HXD/PIN data. 



© 2010 RAS, MNRAS 000, 1-?? 



6 A Lobban et ah 



CM 





^1*1 'l 

ii'i'iUii 


rii,l!{li; 

iiJiltli', 


illtIT- ,viti' ''',4 tiilJl' 

[t'l ^ tlVft 


mm 


'll'HUtli 'llll'- 



0.8 1 1.2 1.4 

Observed Energy (keV) 



Figure 6. A ratio plot showing tlie divergence of tlie XIS 1 detector data 
(liglit grey) witli tliat of tlie XIS 0, 2 and 3 detectors (black) in the soft X- 
ray regime (see Section 3.2). 



To check for consistency, the 22-month time-averaged Swift BAT 
spectra (14-150 keV) were also included (Tueller et al. 2009). This 
provided an extension of the spectrum far beyond the high en- 
ergy limit of the HXD PIN. We restricted the energy range of 
the HXD from 15-50 keV and applied a constant multiplicative 
factor of 1.16 to account for the cross-normalisation at the XIS 
nominal pointing position. The constant in front of the BAT data 
was allowed to be free as the 14-150 keV flux over the 22-month 
period (2004/12/15-2006/10/27) may have varied (time-averaged 
flux Fi4-i50 = 5.31 X 10~^^ ergcm"^ s~^). This provided a 
cross-normalisation factor of 0.75 ± 0.11 for the BAT compared 
to the Suzaku XIS. The power-law component was replaced by 
an exponential cutoff power-law and we note that no cutoff en- 
ergy is required in these data. We constrain the 90% lower limit on 
Ecut > 350 keV. A simultaneous XMM-Newton and BeppoSAX 
PDS observation published by Bianchi et al. (2004) required a high 
energy cutoff with a value of Ecut = 90^20 keV. However, fixing 
the cutoff energy at 90 keV in our Suzaku+BAT spectrum results in 
a worse fit with Ax^ = 35 supportive of the notion that the cutoff 
appears to be at significantly higher energy in these data. 

The residuals in the hard X-ray regime show very little excess 
flux above lOkeV when modelled with a power-law indicating a 
lack of a Compton reflection component. To test for this, we in- 
cluded the PEXRAV model (Magdziarz & Zdziarski 1995) which is 
an additive component incorporating the reflected continuum from 
a neutral slab. We tied the photon index of the power-law incident 
upon the reflector to that of the power-law continuum and fixed 
the elemental abundances to Solar (Anders & Grevesse 1989). We 
also fixed the cosine of the inclination angle of the source to 0.87 
and tied the folding energy to the cutoff energy of the power-law at 
~1000keV, consistent with no cut-off as above. The inclusion of 
the PEXRAV component resulted in a reflection scaling factor value 
of 7? = 0.18lo;i4, consistent with the XMM-Newton/BeppoSAX 
analysis of Bianchi et al. (2004) who find R < 0.19. The inclu- 
sion of this component corresponds to Ax^ ~ 5.0 suggesting that 
this component is only marginally required. This resulted in a fit 
statistic Id.o.f .= 486/446, acceptable with a null hypothesis 
probability of 0.10. 



The next step was to model the entire broad-band spectra by 
including the XIS data below 2.5 keV. As the signal-to-noise ratio 
decreases at lower energies, the spectra were only included down to 
0.6 keV for each XIS. The data were also ignored from 1.7-1.9 keV 
so as not to include the silicon absorption edge due to the detectors. 
The inclusion of these data resulted in a slightly worse fit with a 
null probability of 2.98 x 10^^. Residuals were observed at ener- 
gies < 2keV hinting at the presence of a weak soft excess (e.g., 
Figure 3). Upon closer inspection of these residuals it was noted 
that the XIS 1 (BI) detector data slightly diverged with that of the 
remaining XIS detectors (Figure 6), even when the photon index of 
the power-law continuum was allowed to vary between detectors. 
As the XIS 0, 2 and 3 (FI) were all self-consistent, this divergence 
was possibly caused by calibration effects around the oxygen K de- 
tector edge. To account for this, the data from the XIS 1 detector 
were ignored below 0.7 keV. 

In an attempt to then model the observed soft excess, we added 
a MEKAL thermal plasma component incorporating the emission 
spectra from a hot diffuse gas (Starling et al. 2005). It is worth not- 
ing that a featureless blackbody component models the soft excess 
with an equally good fit as there are no strong lines. However, this 
model was not considered any further as the lack of evidence for 
a strong Big Blue Bump (Wu et al. 1983) suggested little thermal 
emission directly from the disc (Starling et al. 2005). The addition 
of the MEKAL component gave a best-fitting thermal plasma tem- 
perature of k^T = 0.271^5;;;^ keV and resulted in an overall better 
fit with x^ld-o.f- 1104/1022 compared to 1151/1024 before 
the MEKAL component was added"^. The luminosity of the MEKAL 
component was calculated at L = 2.14 x 10'*" erg s^* correspond- 
ing to only 1% of the the total luminosity in the 0.5-10 keV band. 
As no further significant residuals were observed in the spectra 
(Figure 7), this became our accepted broad-band model, the final 
parameters of which are summarised in Table 2 (note that the full 
broad-band model from 0.6-150 keV also resulted in a tighter con- 
straint on the reflection scaling factor, R, the best-fit value of which 
is shown in Table 2). A plot of the relative unfolded model contri- 
butions is shown in Figure 8. 



3.3 Comparison with Past Observations 

3.3.1 XMM-Newton & BeppoSAX 

We tested for any long-term variations in the source by applying 
our best-fit broad-band Suzaku model to the 2001 simultaneous 
XMM-Newton (30 ks exposure) and BeppoSAX PDS (38 ks expo- 
sure) observation (May 2001; see Table 3), as published previously 
by Bianchi et al. (2003, 2004) and Starling et al. (2005). The ra- 
tio of the XMM-Newton EPIC-pn spectrum from 0.3-10 keV to 
the best fit Suzaku model, with continuum parameters described 
in Table 2, is shown in Figure 9. It can be seen that compared to 
the Suzaku XIS spectrum, the pn spectrum is steeper, while overall 
the flux was slightly lower in the XMM-Newton data with a value 
of 2.19 X 10~** ergcm"^ s~* over the 2-10 keV energy range. We 
also note that the flux obtained from the BeppoSAX data was lower 
than that obtained with the Suzaku HXD over the 12-100 keV 



Modelling the soft excess with an APEC component was also consid- 
ered (http://hea-www.harvard.edu/APEC/; Smith et al. 2001) although at 
the Suzaku resolution here, the fit was identical to that obtained with the 
MEKAL thermal plasma. 
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Figure 7. The broad-band spectra of NGC 7213 (in the observed frame) 
with the inclusion of the Swift BAT spectra (blue). A plot of the ratio of the 
residuals to the model desciibed in Section 3.2 and Table 2 is shown in the 
lower panel. No significant residuals are present. 



Model Component 


Fit Parameter 


Value 


1 . Cutoff Power-law" 


r 


1.74 ± 0.01 




High Energy Cutoff 


> 350 




Normalisa tion 




2. Galactic Absorption'' 




1.10 X 10^" 


3. PEXRAV*^ 


\R\ 




4. MEKAL'' 


k^T 


-0.04 




Normalisation 


1 r,Y + 0.59 

-"-■^'-0.38 




Luminosity 


2.14 X 10"'° 


5. Fit Statistics'^ 


xVdof 


1100/1041 




Null Probability 


4.40 X 10-2 



Table 2. Spectral Parameters of continuum fit in the 0.6-150 keV range 
(see Section 3.2 for details). " F, photon index; High energy cutoff in units 
keV; Normalisation in units xlO"^ photons cm—^ s"'^. Local Galactic 
absorption (at z = 0), units cm~^ . \R\, reflection scaling factor, where 
a value K = \ coiresponds to reflection from neutral material subtending 
27rsr. k^T, plasma temperature in units keV; Normalisation (Emission 
Measure = J rienndV) in units xlO^'^cm^''. Reduced chi-squared 
(X^/rfo/) and null hypothesis probability for spectral fit. 




Observed Energy (keV) 

Figure 8. A plot showing the relative contributions of the individual 
model components across the broad-band 0.6-100 keV Suzaiu-l-BAT en- 
ergy range (see Section 3.2 for details). The PEXRAV neutral reflection com- 
ponent is shown in orange over the 1-100 keV energy range whilst the soft 
coUisionally-ionised MEKAL component (at energies <3keV) is shown in 
navy blue. The emission in the Fe K complex is modelled with individual 
Gaussians. The red curve spanning the whole energy range shows the sum 
of all the model components (including the absorbed power-law continuum) 
superimposed on the data. 

range (2.46 x 10"^^ compared to 3.81 x 10^^^ ergcm^^ s"^ re- 
spectively). Since the photon index of the power-law continuum is 
quite hard in this source, this difference could simply arise from 
small changes in the intrinsic power-law. The model is generally in 
good agreement with the data although subtle changes in the con- 
tinuum can be observed, with the spectral curvature being more 
apparent in the XMM-Newton data, e.g., with the spectrum being 
noticeably steeper below 2 keV, but somewhat flatter above 3 keV. 
No strong residuals are present in the iron K band, which suggests 
that the iron line emission has remained constant between the 2001 
XMM-Newton and 2006 Suzaku observations. 



To quantify the changes in the spectrum, the single power-law 
continuum used to fit the Suzaku data in Sections 3.1 and 3.2 was 
replaced with a broken power-law, breaking at 2.19^g'3Q keV with 
r values of 1.84±0.01 and 1.71 ±0.02 below and above this energy 
break respectively. Furthermore, a slight softening of the spectrum 
below 1 keV in the XMM-Newton data, as suggested by the bump 
in the residuals around 0.9 keV (which may due to the Ne IX triplet 
or a blend of emission from iron L-shell lines), indicates that the 
single temperature MEKAL component obtained from the Suzaku 
data was not sufficient to model the soft excess. Thus a second 
MEKAL component was added to the model with a higher tem- 
perature of ksT = OMtoAi 

keV, which significantly improved 
the fit (Ax^ = 19.2) compared to the model with only a single 
temperature plasma. The fit parameters of this best-fit model to the 
XMM-Newton data are summarised in Table 4. 

For consistency this two temperature MEKAL model was then 
applied to the 2006 Suzaku dataset. The temperatures and normal- 
isations of the MEKAL were kept fixed at the best-fit values from 
the XMM-Newton data, as an extended diffuse collisional plasma 
may not be expected to vary significantly over time (note that if the 
parameters are allowed to vary, the values obtained from Suzaku 
are consistent with the XMM-Newton data, within the errors). The 
broken-powerlaw continuum parameterisation was also retained 
from the XMM-Newton fit, although the photon indices and nor- 
malisations were allowed to vary. A comparison of fit parameters 
for the Suzaku and XMM-Newton observations is shown in Table 
4. No other significant variations were observed between the two 
observations; the Fe K line parameters appear to be consistent with 
constant values (within the errors) for the centroid energy and line 
fluxes, while the 2001 BeppoS AX PDS data show no evidence for 
a reflection hump above lOkeV (with R < 0.2; see also Bianchi 
et al. 2004), consistent with what is found by the Suzaku HXD. 
Furthermore, the lower limit on the high energy cutoff value is also 
constrained to >300keV (see also Dadina 2008), in good agree- 
ment with the HXD and Swift BAT. 
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Figure 9. Data/Model residuals of the 2001 XMM-Newton EPIC-pn data 
to the Suzaku model listed in Table 2. The XMM-Newton data show the 
source to have a slightly steeper spectrum compared to Suzaku at energies 
below 2 keV, while overall the continuum flux in the 2-10 keV band is lower 
in the XMM-Newtoji data in 2001 compared to Suzaku in 2006. See Section 
3.3 for details. 



Date 


Mission 


Instrument 


Exposure (ks) 


27/05/2001 


BeppoSAX 


MECS 


61 






PDS 


38 


29/05/2001 


XMM-Newton 


EPIC-PN 


30 


15/12/2004 - 27/10/2006 


Swift 


BAT 


2300 


22/10/2006 


Suzaku 


XIS/HXD 


90 


06/08/2007 


Chandra 


HETG 


148 



Table 3. Log of observations of NGC 7213. See Section 3.3 for details. 



3.3.2 Chandra HETG 

The emission line at 6.60 keV in the Suzaku XIS spectrum is found 
to be rejected at the > 99.9% confidence level as arising from the 
resonance transition of helium-like iron (see Section 3.1 and Fig- 
ure 5). As a consistency check, we modelled the archival Chan- 
dra HETG spectrum (using the latest version of the calibration 
database; v.4.2.2) at FeK with an absorbed power-law and param- 
eterised the emission lines with simple Gaussians consistent with 
the values found by Bianchi et al. (2008) (also see Table 4) with 
emission lines at 6.40, 6.72 and 6.99 keV. Fixing an additional nar- 
row line (a — 10 eV) at 6.60 keV (as required by Suzaku) in the 
HETG spectrum was not required by the data but resulted in a value 
for the equivalent width of EW <21.2eV, which is consistent with 
the equivalent width of EW — 24.4 ± 8.0 eV found in the Suzaku 
XIS. Likewise, including a narrow Gaussian in the Suzaku spec- 
trum with the centroid energy fixed at 6.72 keV (as found by Chan- 
dra) is also not required by the XIS data but yields an upper limit 
on the equivalent width of EW < 14.3 eV, again consistent with the 
equivalent width of EW — 24.0 ± 17.0 eV found by Bianchi et al. 
(2008) with the Chandra HETG. Therefore, it appears that the Fe K 
parameters in both datasets are consistent with each other with no 
evidence of variability detected within the errors. 



4 DISCUSSION 

In this section the possible origin of both the neutral and ionised 
iron K line emission from NGC 7213 is discussed along with its 
implications for the nature of the central engine in this source. 



4.1 The Origin of the Neutral Fe Ka Line 

We first investigated the possibility of whether a distant Compton- 
thick reflector, e.g., such as the pc-scale torus, could account for the 
neutral Fe Ka emission. The tight constraint on the reflection frac- 
tion of i? = 0.09l?);i^, obtained in Section 3.2 (also see Table 2), 
appears to rule out the possibility of the 6.39 keV line originating 
via scattering off Compton-thick matter since, for an FeKa line 
with an equivalent width of ~80eV (as observed here), a strong 
reflection scaling factor value of i? ~ 0.6 would be required 
(George & Fabian 1991). To test this further, the ionised reflec- 
tion model REFLIONX (Ross & Fabian 2005) was used in place 
of the simple 6.39 keV Gaussian emission and PEXRAV model for 
the Compton-scattered continuum off neutral material. The other 
model components, as described in Section 3.2 (also see Tables 
1 and 2), were adopted and are identical in the spectral fits. The 
REFLIONX model consists of the emergent spectrum for a photo- 
ionised optically-thick slab of gas when irradiated by a power-law 
spectrum, with a high energy exponential cut-off of 300 keV, us- 
ing the abundances of Anders & Ebihara (1982). The advantage of 
the REFLIONX model is that it self-consistently computes both the 
reflected continuum and line emission for the astrophysically abun- 
dant elements. 

We initially fixed the iron abundance to Solar, while the red- 
shift of the reflector was found to be consistent with the cosmo- 
logical redshift of the source, with no net (e.g., gravitational) red- 
shift. Given the narrow unresolved iron Kq emission observed in 
the Suzaku spectrum, no additional velocity broadening was ap- 
plied to the reflected spectrum. We also fixed the reflector ion- 
isation parameter'* at a value of ^ = 10 erg cm s~* (the lowest 
value allowed by the model), corresponding to near-neutral iron 
(i.e., iron atoms typically in a low ionisation state corresponding 
to Fe I-XVII). This model provides an upper limit on the reflection 
scaling factor of i? < 0.16, consistent with what was found by 
the PEXRAV model in Section 3, but results in a relatively poor fit 
to the Suzaku+BAT data of x^/d.o.f.= 810/715 (null probabil- 
ity = 7.22 X 10~^). This is due to the fact that the model under- 
predicts the amount of iron Ka emission, leaving a significant pos- 
itive residual at ~6.4keV in the Suzaku XIS data. Allowing the 
Fe abundance to vary to enhance the iron K emission results in an 
acceptable fit of x^ld.o.f.^ 745/714 corresponding to a null hy- 
pothesis probability of 0.21 (with 7? < 0.06). However, in order to 
adequately model the iron Ka line, this requires an overabundance 
of Fe by a factor of ~ 10 with respect to Solar (the 90% confidence 
lower limit on this value is still 4 times Solar). 

It therefore appears that the lack of an observed Compton re- 
flection hump in the data above lOkeV means that the reflection 
and FeKa emission cannot be simultaneously modelled in this 



* Note that in the REFLIONX ionised reflection model the ionisation param- 
eter is defined as f = !kEll and has units erg cm s^^, where F is the illu- 
minating flux incident upon the reflector (integrated over the energy range 
100 eV to 1 MeV) and n is the gas density in cm~^. 
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Model Component 


Fit Parameter 


Suzaku 


XMM-Nemton & BeppoSAX 


Cuanora/HETO 


Flux" 


.Fb 5 2 


1.37 


1.29 






^2-10 


2.44 


2.19 


_ 




^^12-100 


3.81 


2.46 


- 


Broken Powerlaw'' 


Ti 


1 Q4+0-09 

— 06 


1.84±0.01 


- 




-^brcak 


^'■'-'-0.07 


2 1 9+0.34 
^-'^-0.30 


- 




r2 


1.74±0.01 


1.71±0.02 


- 




Normalisation 


6 11+009 


5 92+0.03 
-'■'^—0.04 


- 


MEKAL^ 




FIXED 


0.24+ni 

— U.Uu 






Normi 


FIXED 




_ 




kBT2 


FIXED 




- 




Norm2 


FIXED 




- 




^total 


1.99 


1.99 


_ 


Fe Ka Line*^ 




6.39±0.01 


6.40±0.02 


6.40±0.01 




EWKa 


no ,+11.0 

OJ.l -^Q 7 


85.9lil;« 


120t^° 




Line HuxKa 


lQ+0.28 


ry 00 + 0.38 


9 Q+0.9 
•^■"—0 7 


6.60keVLine<^ 


Eg. 60 


6.6()±0.03 


°-"°-0.05 


g 72+0:01 

— 0.(32 




EWe.eo 


24.1±7.9 


27.4±11.6 


24± 17 




Line Huxg.eo 


0.67±0.22 


0.78±0.33 


0.7 ±0.5 


6.95 keV Line'' 


£^6.96 


6.95±0.03 


6.91±0.08 


-i_n no 

f. QQ+O.U^ 




EW6.96 


62 3+1®° 


37.1±18.8 


60 ±30 




Line HUX6.95 


1 40+0-36 


0.89±0.45 


1.3 ±0.6 


Statistics 


X^fd.o.f. 


747/723 


846/872 






Null Probability 


0.26 


0.73 





Table 4. Spectral parameters of the broad-band fit in the 0.6-1 50 keV energy range for the 2006 Suzaku and 2001 simultaneous XMM-Newton and BeppoSAX 
observations (see Section 3.3 for details). The Fe-hne complex parameters from the 2007 Qiandra/HETG observation (reported by Bianchi et al. 2008) are also 
shown for ease of comparison. " Continuum flux in the specified range in units xlO"^^ ergcm"^ s~i. F, photon index; Ebreak in units keV; NormaUsation 
in units x 10~3 photons cm~^ s"^. k^T, plasma temperature in units keV; NormaUsation (= J n^riYidy) in units x 10®^ cm~3; Total mekal luminosity 
from 0.5-lOkeV in units 10*0 ergs"^ . AH MEKAL fit parameters in the Suzaku dataset are fixed at the best-fit values from the XMM data. Bune. line 
energy in units keV; EW, equivalent width in units eV; Line Flux in units xlO~^ photons cm"^ s"^. 



way, seemingly ruling out a reflection origin for the Fe Ka emis- 
sion, as also suggested by Bianchi et al. (2003) on the basis of the 
BeppoSAX data. Indeed, an acceptable fit can only be obtained 
if the abundances are assumed to be greatly super-Solar, at odds 
with the modest spread of values found by Perola et al. (2002) 
from a sample of nine bright Type 1 Seyferts and NELGs (Narrow 
Emission-Line Galaxies) observed with BeppoSAX. Thus is ap- 
pears unlikely that the 6.39 keV emission originates via reflection 
off Compton-thick matter unless the material covers a very small 
solid angle (< 1 sr) and is extremely iron over-abimdant. 

Instead it is perhaps more likely that the near-neutral iron Ka 
line originates in Compton-thin matter, covering a higher fraction 
of 47r steradians solid angle. Indeed, an estimate of the column den- 
sity of the Ka-emitting material can be made using the calculations 
of Yaqoob et al. (2010) where an analytic expression relating the ef- 
ficiency of the Fe Ka line production and the column density of the 
emitting material is derived in the optically-thin limit. The produc- 
tion efficiency of the Fe Ka line is calculated by: 

XFeK. - .00 ^_p^^- (1) 
JfcK 

Here, Ek is the threshold energy for FeK-shell absorption 
and /FeKa.n refers to the line flux renormalised to an incident 
continuimi with a flux of 1 photon cm~^ s~ikeV~i at IkeV. T 
is the photon index assuming an incident power-law continuimi. 
Adopting the Vemer et al. (1996) value for of 7.124keV, we 
calculate an FeKa line production efficiency of ~1% for NGC 



7213. In the Compton-thin case, we find that this results in an 
estimate on the column density of the Ka-emitting material of 
A^H 2 X 10^^ cm"'^ using the analytic expression derived by 
Yaqoob et al. (2010) (equation 4 in the aforementioned publica- 
tion). Although this value is consistent with that found by Bianchi 
et al. (2008), the expression is valid only in the Compton-thin limit 
which begins to break down for iVn > 2 x 10^^ cm~^ as the op- 
tical depth of the FeK line photons to absorption and scattering 
becomes non-negligible (see Yaqoob et al. 2010, Figure 2). How- 
ever, accounting for these effects, according to the calculations of 
Yaqoob et al. (2010) and Murphy & Yaqoob (2009), a column den- 
sity of A'h ~ 3 — 4 X 10^^ cm~^ can result in an FeKa line effi- 
ciency of 1 % for a face-on geometry covering 2n sr in their toroidal 
X-ray reprocessor model. Thus the Fe Ka line may originate in a 
Compton-thin torus or perhaps the outer BLR clouds, as suggested 
by Bianchi et al. (2005), although the covering fraction would per- 
haps be sUghtly high in the latter case (Netzer & Laor 1993). 

4.2 The Origin of the Highly Ionised Fe 

We investigated the possibility of an ionised accretion disc as a 
potential origin for the 6.60 and 6.95 keV emission lines. To test 
this scenario, we attempted to model the two ionised emission lines 
by a REH.IONX component with a high ionisation parameter of 
log^ = S.OIq i ergcms"'^. The remainder of the Suzaku spec- 
trum was modeUed as before, i.e., a cut-off power-law for the con- 
tinuum emission, a single temperature thermal mekal component 
for the weak soft X-ray excess and a narrow Gaussian centered at 
6.39 keV for the neutral iron Ka emission. No relativistic blurring 
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Figure 10. A plot showing the inability to model the highly ionised Fe emis- 
sion at 6.60 keV and 6.95 keV with an ionised reflector therefore suggesting 
a non-reflection origin for these lines (as described in Section 4.2). The con- 
tinuum level is shown in green and the lower panel shows the ratio of the 
residuals to the model. The neutral Ko line at 6.39 keV is modelled with a 
Gaussian. The vertical dotted lines show the expected line energies of, from 
left to right, Fe K«, Fe XXV forbidden, Fe XXV ls-2p resonance, Fe XXVI 
ls-2p and Fe K/3 in the observed frame. 




6 6.5 7 



Observed Energy (keV) 

Figure 11. A plot showing the fit obtained when the highly ionised Fe is 
modelled with a photo-ionised XSTAR grid with a = 3000 kms~^ and 
log^ = 4.4ergcms~^ (as described in Section 4.2.2). The continuum level 
is shown in green and the lower panel shows the ratio of the residuals to the 
model. The neutral Fe Ka line at 6.39 keV is modelled with a Gaussian. The 
vertical dotted lines show the expected line energies of, from left to right, 
FeKa, FeXXV forbidden, FeXXV ls-2p resonance, FeXXVI ls-2p and 
Fe K/3 in the observed frame. 

was applied. This resulted in a poor fit of Id.o.f .= 791/717 
(null probability = 2.79 x 10^^) with the model unable to account 
for the ionised emission from either Fe XXV or Fe XXVI (see Fig- 
ure 10). Indeed due to the large intrinsic electron scattering depth 
at such high values of ^, the lines in the reflection model become 
too broadened and so are unable to model the relatively narrow 
6.60 and 6.95 keV lines that are observed in the Suzaku spectrum. 
Furthermore the 90% confidence upper-limit on the reflection frac- 



tion for the highly ionised reflector is restricted to 7i < 0.06, 
which seemingly allows us to reject the presence of a highly ionised 
Compton-thick medium in NGC 7213, such as a highly ionised, but 
Compton-thick inner accretion disc. 

In addition, the inability to model the neutral Fe Ka emission 
with a low ionisation reflector (see Section 4.1) coupled with the 
apparent lack of a Compton hump > 10 keV means that this source 
exhibits no evidence at all for any Compton-thick material, either 
neutral or ionised. Furthermore, the data do not appear to show any 
relativistic signatures since the neutral and ionised iron lines appear 
to be narrow (or only moderately broad in the case of Fe XXVI), 
thus ruling out any emission from the innermost regions around 
the black hole. Combining the lack of evidence for any Compton- 
thick matter in this source with the lack of any relativistic signature 
seems to suggest the complete absence of an inner optically-thick 
accretion disc in NGC 7213. Instead it may support the notion of 
an accretion disc which is truncated at some radius with the inner 
regions replaced by some form of Radiatively Inefficient Accretion 
Flow (RIAF; Narayan & Yi 1995), where the low efficiency of the 
inner, hot corona leads to much of the energy being advected into 
the black hole rather than radiated away. 

4.2. 1 CoHisionally-Ionised Model 

We instead investigated the possibility of a coUisionally-ionised 
origin for the highly ionised FeK emission which we modelled 
with a high temperature broadened APEC component (Smith et al. 
2001). This required a temperature of fceT = 13.03l2'g6 ^^"^ 
velocity broadening, a ~ 2500 km s~^, resulting in a fit statistic 
of x^/d.o.f.— 788/716. Inspection of the data around the iron K 
band reveals that the collisionally-ionised model over-predicts the 
contribution of the resonance line compared to the forbidden line of 
Fe XXV, which leaves some excess line flux near 6.60 keV unmod- 
elled, although the FeXXVI line is modelled well. Thus, although 
the collisional model cannot be ruled out with high confidence in 
the present data, given that the energy of the 6.60keV line appears 
more consistent with the forbidden rather than resonance transition 
of Fe XXV, an additional photo-ionised component would perhaps 
be required to model this line (see Section 4.2.2). 

An estimate on the density of the collisionally-ionised mate- 
rial can be calculated from the normalisation of the APEC code, 
which returns J dV ~ lO^^cm"^. If the observed broadening 
of the line is assumed to be due to a Keplerian orbit (ct ~ 0.1 keV 
corresponding to a few 10'* km s~* ; see Table 1), we can calculate 
an estimate on the radius of the emitting material from the central 
black hole of a few 10'' Rg, where 1 Rg = 1.5 x 10*^ cm for a black 
hole mass of 10® Mq (Nelson & Whittle 1995). Adopting this ra- 
dius of i? ~ 1.5 X lO^'^ cm results in an estimate of n ~ 10^ cm~^ 
for the density of the emitting material. 

If the collisional gas is indeed responsible for some of the 
highly ionised Fe emission then we note that the conditions of 
the gas must be different from those in the lower temperature gas 
used to model the weak soft excess at energies < 1 keV (see Sec- 
tion 3.2). The two components require very different temperatures 
(0.271|5;m and 13.03l2;86 keV for the low and high temperature 
gases respectively) and we note that any contribution from one dis- 
tinct zone of emission to the other is negligible (i.e., no significant 
emission from the lower temperature zone is found to be contribut- 
ing the emission at Fe K and vice-versa). However, the two zones 
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of gas could be linked by their origin albeit on completely different 
scales with the emission from the lower temperature gas originating 
at a much greater distance from the black hole. 



4.2.2 A Photo-ionisation Model for the Ionised Iron K Emission 

An alternative origin for the lines could be emission from a photo- 
ionised, but Compton-thin plasma. We modelled this scenario using 
the XSTAR 2.11nll code (Kallman & McCray 1982), which incor- 
porates the abundances of Grevesse, Noels & Sauval (1996). We 
initially modelled the lines using a single zone of emission with So- 
lar abundances and a turbulent velocity width of 200kms~^. The 
best-fit value of the ionisation parameter was log ^ = 3.7 ± 0.3 
(where the ionisation parameter in xstar is defined as in equation 
3), which resulted in an acceptable fit with Id.o.f .— 779/715 
(null probability = 4.73 x 10~^). Even so, some slight excess resid- 
uals were apparent, particularly around the 6.95 keV line, which 
might suggest there is some intrinsic velocity broadening of the 
ionised lines. For instance, as measured from Section 3.1, the H- 
like line appears to have a IWHM of ~ 1 000 km s ~ ^ . 

Thus an alternative xstar grid was used with a higher turbu- 
lence velocity of a = 3000 km s~^. This resulted in a slightly bet- 
ter fit of x^!d.o.f.= 765/715 (null probability = 8.50 x 10"^), 
with an ionisation parameter of log 5 = 4.4±0.1 (see Section 4.3.3 
for a calculation of the lower limit on the column density). This fit 
is significantly better than the one obtained with the coUisionally- 
ionised model (Section 4.2.1) with Ax^ ~ 23 between the two 
models. A plot of this model superimposed on the data is shown 
in Figure 11. We also calculated the 90% uncertainty on the red- 
shift of the zone allowing us to constrain the net velocity shift of 
the ionised emitter to be v = -t-650li5oo km s~^, where a positive 
velocity denotes rcdshifted / infalling material. Thus the data do 
not formally require a velocity shift in this model. Further zones of 
ionised matter are also not required by the data. 

Both high and low turbulence models appear to give good fits 
over the iron K band, with the 6.95 keV emission line originating 
from the Fe XXVI Lya doublet and the 6.60 keV He-like line aris- 
ing due to a blend of the forbidden and intercombination lines at 
6.637 and 6.668 keV respectively. In the low density limit assumed 
in the XSTAR model here (where n < 10^® cm~^), the He-like 
emission is from an approximately equal mixture of the forbid- 
den and intercombination lines, with a negligible contribution from 
the resonance line at 6.700 keV. Note that some weak emission via 
satellite lines of lower ionisation iron (i.e., FeXXIII-XXIV) could 
also be contributing to this blend of emission although at the high 
ionisation parameter inferred here (of log ^ ^ 4), this contribu- 
tion is likely to be negligible (Kallman & McCray 1982). However, 
future calorimeter-based spectroscopy, e.g., with Astro-H, will be 
required to spectrally resolve all the line emission components as- 
sociated with the He-like triplet of iron and to constrain the intrinsic 
velocity broadening. 

We also note that if the highly ionised iron lines are originat- 
ing in photo-ionised gas then the possibility remains that a simi- 
lar photo-ionised gas but with a much lower ionisation parameter 
could be responsible for the weak soft X-ray lines observed at ener- 
gies < 1 keV. However, Starling et al. (2005) analysed the results of 
an XMM-JVewton RGS observation of NGC 7213 and found that a 
coUisionally-ionised thermal plasma was preferred by the data from 



consideration of the ' G' ratio (Porquet & Dubau 2000) of the O Vii 
triplet. 

4.2.3 The Location of the Highly Ionised Gas 

From consideration of various photo-ionised and collisionally- 
ionised models, it appears that the 6.60 and 6.95 keV lines are con- 
sistent with originating in a photo-ionised medium, although the H- 
like Lya fine appears to be quite strong with an equivalent width of 
~60 eV. Bianchi & Matt (2002) calculate the equivalent widths of 
Fe XXV and Fe XXVI lines with respect to both the reflected and to- 
tal continua and show that an equivalent width of ~15 eV for H-like 
Fe would be expected if it originated in material of column density, 
TVh ~ 10^^ cm~^, with an ionisation parameter of log^ ~ 3.5, 
photon index of F = 1.7 and a covering fraction, / = 1. To test 
whether our cquivalcnl width of ~60 eV was feasible we calculated 
a lower limit on the column density of the photo-ionised material 
in the case of one single zone of emission. 

For a uniform, spherical, Compton-thin shell, the normalisa- 
tion of the photo-ionised emission component is defined within the 
XSTAR codeby^: 

k ~ fcov (2) 

where fcov is the covering fraction of the material (i.e., fcov = 
1 for matter covering An sr^^), Lion is the ionising luminosity in 
units 10"*** erg s^'^ from 1 to 1 000 Rydbergs and D is the distance 
to the source in kpc. Physically the XSTAR normalisation is simply 
proportional to the observed X-ray flux of the source multiplied 
by the covering fraction of the photo-ionised gas. Thus the appro- 
priate value of this normalisation, k, for NGC 7213 can be cal- 
culated assuming a covering fraction of 1. For the luminosity and 
distance of NGC 7213, this results in a value of A; = 9.0 x 10"^ 
in units 10^* ergs~^ kpc"'^, where we have adopted a distance of 
25 000 kpc to NGC 7213 (Mould et al. 2000). The value used for 
the ionising luminosity was derived from an extrapolation of the 
broken power-law continuum in Table 4, integrated from 1 to 1 000 
Rydbergs and is found to be 5.61 x 10*^ ergs~^, comparable to 
the estimate of the bolometric luminosity of ~9 x 10*^ ergs~^ by 
Starling et al. (2005). 

In the model used in Section 4.2.2, we have fixed the nor- 
malisation of the additive xstar emission component to the above 
value and thus assumed a fully covering shell of gas around 
NGC 7213, with the appropriate source luminosity and distance 
from above. The column density was then allowed to vary within 
the xstar model, in order to fit the ionised iron K emission lines. 
This provided a best fit value for the column density of A^h = 
4.0lo;g ^ lO^'"* cm"^. The 90% confidence lower Umit on the col- 
umn density is A^h > 3.2 x 10^"* cm~^ for a fully covering spheri- 
cal shell of gas, with an ionisation parameter of log ^ ~ 4 as above. 
For a shell that does not fully cover the source, then the column 
density will need to be higher to compensate for the lower cov- 
ering. This is in good agreement with the calculations of Bianchi 
& Matt (2002), who predict an equivalent width for Fe xxvi Lya 
of 15 eV for a column density of 10^^ cm~^, compared to the ob- 
served 60 eV equivalent width in the case of NGC 7213, but for a 

^ see 

http://heasarc.gsfc.nasa.gov/docs/software/xstar/docs/htinl/node94.html 
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column density approximately 3^ times higher. 

In order to better constrain the origin of the Fe XXVI Lya line 
it is important to estimate a value for the distance of the emitting 
material. Assuming a uniformly ionised, spherical shell of gas, the 
ionisation parameter in xstar is defined as: 



(3) 



and has units erg cm s ~ ^ where Lion is the ionising luminosity 
from 1 to 1 000 Rydbergs, n is the gas density in cm~^ and R is the 
radius of the absorbing / emitting material from the central source 
of X-rays. Combining this with the column density which is given 
by: 

poo 

Nh= ndR (4) 
yields an estimate on the inner radius of the emitting material: 



(5) 



5 X 10*^ erg s ^ (derived above). 



3 X lO'^^ cm then gives a lower limit on R 



Assuming values of Lj 
C ~ 5 000 and iVn 
of ~3 X 10^'' cm corresponding to a value of ~200 Rg (assuming 
a black hole mass of 10*^ Mq; Nelson & Whittle 1995) and infers 
an electron density of rie ~ 10® cm~^ (equation 3). This radius 
is also consistent with the FWHM of 10 000 km s"^ of the H- 
like Fe line which, if the broadening is assumed to be the intrinsic 
broadening due to a Keplerian orbit, provides an estimate on the 
Une emitting radius of a few 10^ Rg. 



4.3 NGC 7213 as a Low Luminosity AGN 

4.3.1 The Inner Advective Flow in NGC 7213 

From 2-10 keV, the X-ray spectrum of NGC 7213 resembles that 
of a typical Type 1 Seyfert Galaxy where the spectrum is domi- 
nated by a power-law continuum of F = 1.75 and near- neutral 
FeKa emission at 6.39 keV. This neutral Ka emission may orig- 
inate from Compton-thin material of A^h ~ 3 — 4 x 10^^ cm~^ 
possibly in the outer BLR or a Compton-thin torus (Bianchi et al. 
2008). We also detect significant emission from highly ionised ma- 
terial located close to the central source with Fe xxv and Fe xxvi 
perhaps originating in a photo-ionised medium with a column den- 
sity A^H ~ 3 X 10^-^ cm invoked to match the high observed 
EW of ~60eV. This emission is likely to originate at a distance of 
J? ~ 10^ - 10* Rg from the black hole / X-ray source. 

Given the tight constraints on reflection from both neutral and 
ionised material {R < 0.16 and R < 0.06 respectively), the lack of 
any relativistic signatures and the very weak Big Blue Bump (Wu 
et al. 1983; often interpreted as thermal emission from the disc), 
this appears to rule out the presence of an inner, 'classic' optically 
thick, geometrically thin accretion disc (Shakura & Sunyaev 1973) 
envisaged in the unification scheme of AGN (Antonucci 1993). 
Instead, we suggest that the accretion disc is maybe truncated at 
some radius on the order of 10^ — 10* Rg (~ 0.01 pc) with the in- 
ner regions perhaps replaced by a Radiatively Inefficient Accretion 
Flow (RIAF; Narayan & Yi 1995) consisting of highly ionised, low 



density (n ~ 10® cm~^), Compton-thin gas covering some signifi- 
cant fraction of 47r sr. In this scenario, the low accretion rate of the 
source (0.07% L-^dd), perhaps due to a lack of available accreting 
material, may not allow the infalling material to cool sufficiently 
in order for a standard thin accretion disc to form. Instead, a stable 
accretion flow can still occur if the material takes on the form of an 
optically-thin, hot corona, where most of the material is advected 
across the event horizon as opposed to radiating away the energy it 
has acquired in moving close to the black hole. 

Such a hot, inner flow is expected to comprise of a low- 
density plasma whereby coUisional processes dominate over photo- 
ionisation due to the high temperature. At radii below about 
lOORg, the ion and electron temperatures diverge forming a two- 
temperature medium with Te ~ 10® — 10^° K and Tion approaching 
lO^'^ K in the innermost regions (Narayan & Yi 1995). The elec- 
tron temperature To is then expected to fall as 10^^ (K) /R for R 
> 10^, where R is in Schwarzschild units. However, further out 
at R ~ 10^ - 10" Rs, a plasma temperature of a few 10* K (i.e., 
KbT ~ ISkeV, as observed in the coUisional APEC model), cor- 
responding to i? ^ 10* Rs, would produce emission from both 
He-like and H-like Fe, as observed. However, while the Fexxvi 
emission could plausibly originate in such a plasma, the presence 
of the Fe XXV forbidden line suggests that the He-hke Fe cannot be 
solely produced in such a collisionally-ionised medium. Thus an al- 
ternative picture could instead be that both the Fc XXV and Fe XXVI 
emission lines originate in photo-ionised gas, perhaps at the transi- 
tion region between the RIAF and the cold, outer accretion disc at 
a radius R ~ 10^ — 10* Rg from the black hole. Alternatively, the 
Fe XXV / XXVI emission lines may arise from a hybrid of photo- 
and collisionally-ionised processes. 



4.3.2 NGC 7213 as a Low/Hard State Source 

The accretion rate of NGC 7213 is much lower than the pre- 
dicted "critical" value of -^2% LEdd (Maccarone 2003) whereby 
the high/soft state in X-ray binaries can be observed. As a result, 
one interesting possibility is that NGC 7213 is an AGN analogue 
of the low/hard state observed in Galactic Black Hole Candidates 
(GBHCs). Long-term monitoring of NGC 7213 with RXTE (Phil 
Uttley, priv. comm.) shows that the AGN is only slowly variable, 
indicating a relatively low frequency PSD break. This would be 
consistent with NGC 7213 having a relatively high black hole mass 
(e.g., Mbh = 10*^ Mco; Nelson & Whittle 1995; Woo & Urry 2002) 
but a low accretion rate compared to Eddington, consistent with 
the scaling relations in the timing properties seen between AGN 
and GBHCs (McHardy et al. 2006). This is further supported by 
the SED of NGC 7213 which suggests that this object has inter- 
esting radio properties lying between those of Radio-Loud and 
Radio-Quiet Quasars. Indeed, taking the 5 GHz radio flux and B 
band flux (host galaxy-subtracted) of Sadler (1984) and Halpem 
& Filippenko (1984) respectively and using the equation for radio- 
loudness^' of Wilkes & Elvis (1986) gives a value of _Rl ~ 2 sug- 
gesting that NGC 7213 is intermediate between Radio-Quiet AGN 
and Radio-Loud AGN such as Radio-Galaxies and Blazars. How- 
ever, Panessa et al. (2007) also find that a radio-loudness of ~2 is 



® i?L = logioC^S GHz/^b), where Fr, caIz and Fb are the 5 GHz and B 
band fluxes respectively. Typically, a value of i?L > 1 signifies a Radio- 
Loud object. 
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not so uncommon in Seyfert galaxies. 

The inability to constrain the high energy cutoff of the X-ray 
emission (i.e., >350keV) could indicate that the continuum emis- 
sion may be of non-thermal origin with one possibility being that 
some of the hard X-ray emission that we observe with Suzaku orig- 
inates from the base of a jet'^. Hameed et al. (2001) imaged NGC 
7213 in the optical band and discovered a giant Hq filament ap- 
proximately 19kpc from the nucleus. They suggest that such a fil- 
ament could be the signature of neutral gas shock-ionised by the 
interactions of a jet. A more recent 8.4 GHz Long Baseline Ar- 
ray (LBA) radio study of NGC 7213 (Blank et al. 2005) reports 
that the source is unresolved on the scale of 3 milli-arcseconds 
(corresponding to ^10** Rs at the distance and black hole mass of 
NGC 7213), just showing a core, indicating that the jet could be 
orientated face-on. At lower frequencies, there is also evidence for 
a large-scale structure (30-40"; Blank, Harnett & Jones 2005 and 
reference therein) which could possibly be a signature of the ex- 
tended radio lobes. Consequently, NGC 7213 is perhaps consis- 
tent with the hypothesis of Falcke et al. (1996) whereby the radio- 
intermediate objects are similar to Radio-Quiet AGN but with mod- 
erate beaming from pc-scale jets orientated face-on to the observer. 



4.3.3 The Origin of the High Energy Continuum 

NGC 7213 is not detected to date with the Fermi LAT gamma-ray 
instrument at ~ GeV energies (Abdo et al. 2010) where the inverse- 
Compton emission from the jet would be expected to dominate. 
The Fermi LAT flux limit corresponding to the detection threshold 
of Abdo et al. (2010) at the Galactic co-ordinates of NGC 7213 and 
forT = 1.75 is Fo.i-iooGeV < 3 x 10"^ photons cm"^s"\ Ex- 
trapolating our best-fit broad-band Suzaku model to GeV energies 
over-predicts the 7-ray flux by a factor of ~ x 100 returning a pre- 
dicted photon flux of Fo. 1-100 GeV ~ 3 x 10~^ photons cm^'^ s^^. 
This implies that the X-ray continuum does in fact roll over at en- 
ergies >350keV. 

In order to be consistent with the 0.1-100 GeV flux limit from 
Fermi, we require that the E-folding energy of the power-law com- 
ponent must be i?cut < 25 MeV. This, combined with the lower 
limit on the high-energy cutoff from the combined Suzaku I Swift 
data means that 350 keV <Ecut <25MeV, consistent with the 
predicted electron temperature of the hot, inner flow (see Section 
4.3. 1). This suggests that thermal Comptonisation is responsible for 
the X-ray continuum and that any non-thermal contribution from 
the inverse-Compton component associated with a jet may be neg- 
ligible in this source. 

Furthermore, the high EW of the observed emission lines may 
also suggest that there is very Httle dilution of the X-ray contin- 
uum by a jet. We note that other radio-loud sources such as the 
Broad-Line Radio Galaxies (BLRG) 3C 120 and 3C 390.3 do also 
show fairly strong Fe K line emission with EWs on the order of 50- 
100 eV. However, in the case of 3C 390.3, Sambruna et al. (2009) 
argue from the Suzaku data and the overall radio - 7-ray SED that 



' Note that the high energy cutoff was measured with BeppoSAX in 2001 
to be i?cut = SOtjo keV (Bianchi et al. 2004). However, we find no evi- 
dence for a cutoff energy <350keV with our combined Suzaku XIS, HXD 
and Swift BAT spectra. 



the jet makes a minimal contribution to the X-ray continuum emis- 
sion. Likewise, Kataoka et al. (2007) argue a similar case from the 
Suzaku observation of 3C 120 and conclude that the putative jet 
component does not dilute the FeK emission. In addition, com- 
paring the ratio of the 1-100 GeV 7-ray flux from Fermi to the 
2-10 keV X-ray flux for NGC 7213 with that of 3C 111 (the only 
BLRG detected by Fermi to date) and 3C 273, we find that the ra- 
tio is higher by a factor of > x 6 for the two radio-loud AGN. By 
comparison, like for NGC 7213, none of the X-ray bright Type 1 
Sej^erts appear to have been detected with Fermi thus far. 

As a further test, we did attempt to model the X-ray contin- 
uum with a double power-law component consisting of a hard spec- 
trum to model any possible emission from a jet and a much softer, 
Se3^ert-like spectrum to model the photo-ionising nuclear X-ray 
emission. Statistically speaking, this fit is not required by the data 
as it only yields an improvement of Ax^ ~ 2 for an additional two 
parameters of interest. Upon fixing the photon index of the softer 
power-law at F = 2, we find that the normalisation of this compo- 
nent becomes very small with an upper-limit corresponding to just 
10% of the normalisation of the main power-law. So it seems that 
a double power-law model is not required by the data and that the 
X-ray continuum is best represented by a single power-law compo- 
nent. 

The photon index of the X-ray continuum has a best-fit value 
ofF = 1.75±0. 02 and therefore has only a slightly flatter spectrum 
than the tj'pical values usually associated with Radio-Quiet Quasars 
(RQQs) and Type 1 Seyferts. For example. Reeves & Turner (2000) 
find a mean value of F = 1.89 ± 0.05 from a sample of 27 RQQs 
observed with ASCA (Advanced Satellite for Cosmology and As- 
trophysics) and Nandra & Pounds (1994) find a mean value of 
F = 1.95 with a dispersion of cr = 0.15 from their sample of 
Seyfert galaxies observed with Ginga. Porquet et al. (2004) also 
find a mean value of F = 1.90 with a dispersion of a = 0.27 from 
a sample of 14 RQQs observed with XMM-Newton. The photon in- 
dex of NGC 7213 is, however, consistent with those found in other 
low-luminosity AGN such as M81 (Young et al. 2007) and NGC 
4579 (Dewangan et al. 2004). Interestingly, the spectrum in NGC 
7213 does appear to be somewhat steeper than the predicted photon 
index of F = 1.4 from the relation between mass accretion rate and 
photon index of Papadakis et al. (2009) given its calculated accre- 
tion rate of 0.07% Ledd (an accretion rate of ~2% Ledd would 
be required to obtain F = 1.75). So it seems that NGC 7213 may 
not strictly follow the positive correlation between spectral steep- 
ness and accretion rate in AGN and X-ray Binaries suggested by 
Shemmer et al. (2006) although the high energy cutoff of the X- 
ray continuum (i.e., 350 keV <i?cut <25 MeV) observed here may 
suggest that the hard X-ray emission that we observe with Suzaku 
and Swift could be consistent with a very hot inner-flow, compati- 
ble with NGC 7213 having a low mass accretion rate. 



5 CONCLUSIONS 

NGC 7213 is an unusual AGN as it consistently exhibits no evi- 
dence for a Compton reflection component unlike in other Type 1 
Seyferts (Perola et al. 2002; Dadina 2008). The time-averaged con- 
tinuum emission is well fitted by a single power-law of F = 1.75 
and from consideration of the combined Suzaku and Swift BAT 
data and the Fermi flux limit we constrain the high energy cutoff to 
be 350 keV <Ecut <25 MeV. The limits on reflection for the neu- 
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tral and ionised cases from the reflionx model we. R < 0.16 and 
R < 0.06 respectively suggesting that a significant Compton-thick 
reflector (e.g., from the inner disc or Compton-thick torus) is ab- 
sent in this source, consistent with previous findings (e.g., Bianchi 
et al. 2003). Nonetheless, a significant FeK complex is observed 
above 6 keV appearing only in emission. The line from neutral Ka 
dominates (6.39 keV; EW ^ 80 eV) with further contributions from 
FeXXV and FeXXVlLya (6.60 and 6.95 keV respectively; also 
see Starling et al. 2005, Bianchi et al. 2008). Furthermore, in this 
observation we also find that the FexxviLya emission appears 
to be somewhat resolved in the Suzaku spectrum with a FWHM 
~ 10 000 km s^^ and that the emission from FeXXV appears to 
be consistent with the forbidden transition from helium-like iron at 
~6.64 keV as opposed to the resonance transition at ~6.70 keV. 

The neutral Ka emission cannot be modelled via reflection 
off Compton-thick matter. However, an origin in a Compton-thin 
plasma covering a significant fraction of An sr is feasible with an in- 
ferred column density of A'h ~ 3— 4x 10^^ cm~'^, again consistent 
with the findings of previous observations with Chandra and XMM- 
Newton. Likewise here, the emission from highly ionised iron can 
also be modelled with a substantial colurrm (A'h ~ 3 x 10^^ cm~^) 
of photo-ionised matter if a location close to the central engine is 
invoked to explain the inherent broadening and the high ionisation 
state. Given the lack of either neutral or ionised reflection coupled 
with the apparent absence of any relativistic signature in the spec- 
trum, it appears that an inner, optically-thick accretion disc may be 
absent in this source. Instead, the accretion disc in NGC 7213 is 
most likely truncated at some radius with the inner regions perhaps 
replaced by an advective accretion flow (e.g., RIAF; Narayan & Yi 
1995). The FeXXV / XXVI emission could then be the ionised sig- 
nature of such a hot, optically-thin plasma originating in material a 
few 10^ Rg from the central X-ray source. 
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